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Section  A.  Project  Summary 

This  is  the  SEAKEY  final  report  for  the  BBN-RVS  team  for  the  Phase  I  reporting  period, 
February  2014  through  February  2015.  Our  tasks  under  this  project  include  identifying  key 
challenges  to  transition  QKD  to  a  naval  environment  and  determination  of  optimum  operating 
parameters  for  a  naval  QKD  system,  in  pursuit  of  the  Seakey  end-of-program  goal  of  achieving 
key  rates  of  1  Mbit/sec  at  a  range  of  10-30  km  in  a  naval  environment. 

Program  Highlights: 

1 .  Fundamental  upper  bound  to  the  key  rate  achievable  using  any  QKD  protocol  over  a 
lossy  channel. 

2.  Developed  a  MATLAB  tool  for  the  evaluation  of  secret  key  rates  under  turbulent 
propagation  (a  theoretical  model)  and  atmospheric  absorption  and  scattering  (a 
numerical  model  taken  from  MODTRAN),  while  employing  a  single  spatial  mode. 

3.  A  quantitative  trade  study  of  using  multiple  spatial  modes  and  finding  up  to  how  many 
spatial  modes  may  yield  a  perceivable  key  rate  benefit. 

4.  Development  of  a  quantum-secure  direct  communication  protocol  with  near-optimal 
rate-loss  tradeoff  that  uses  laser  light  modulation  and  homodyne  detection,  that  is 
immune  to  a  quantum-limited  passive  eavesdropper. 

5.  Identification  of  a  potentially  far-reduced-complexity  high-rate  CV  KD  protocol. 


This  report  is  organized  as  follows: 

1.  This  document  (containing  minor  programmatic  points,  organizational  notes  and 
summary  of  progress  in  Phase  I); 

2.  A  PPTX  document  (SEAKEY_RESULTS_Y1)  containing  the  summary  presentation 
of  Phase  I;  and 

3.  A  folder  (Model)  containing  the  GUI  m-file  lnput_parameters.m  used  to  calculate  key 
rates  with  inputs  from  the  MODTRAN  database. 


Section  B.  Technical  Progress 

In  this  section,  we  describe  the  Statement  of  Work  (SoW)  proposed  for  Phase  I  and  our 
progress  against  each  of  the  tasks  in  the  SoW. 

Task  1  Identify  key  challenges  for  transition  QKD  to  a  naval  environment.  Examine  rate- 
distance  tradeoffs  with  turbulence,  scattering,  absorption,  background  etc. 
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SoW:  The  first  step  addressing  this  problem  will  be  to  identify  all  of  the  vulnerabilities  QKD  will 
be  subject  to  when  deployed  in  a  naval  environment.  The  evolution  toward  designing  a  robust, 
high-performance  QKD  system  for  a  naval  environment  will  begin  with  a  study  to  identify  the 
key  challenges  to  QKD  in  such  an  environment,  and  quantitatively  assess  their  impact  on 
system  performance.  This  initial  focused  four  month  analysis  will  categorize  the  effects  most 
deleterious  to  fielding  a  high  performance  QKD  system. 

Phase  I  Results: 

The  Seakey  team  has  focused  on  investigating  the  rate-loss  behavior  of  the  two  classes  of 
quantum  key  distribution  protocols  -  discrete-variable  (DV)  and  continuous-variable  (CV).  In 
recent  work  (arXiv:1310.0129),  Guha  et  al  defined  the  Rate-Upper  Bound,  which  sets  an  upper 
limit  on  the  secret  key  capacity  of  a  pure-loss  bosonic  channel.  As  seen  from  Figure  1,  CV  and 
DV  protocols  have  the  same  optimal  rate-loss  scaling,  given  by  R~q.  In  particular,  CV  binary- 
phase-shift-keying  (BPSK)  and  DV  (polarization,  or  time-bin  encoded)  BB84  without  decoy 
states,  both  yield  a  worse  (R~r|2)  scaling.  Thus,  an  extension  of  the  CV  BPSK  protocol  with  a 
few  additional  modulation  levels  (but  far  fewer  from  a  QAM-sampled  discretization  of  the  full 
Gaussian  distribution  of  amplitude  and  phase,  that  CV  demonstrations  use)  should  retrieve  the 
optimal  (R~q)  key  rate  scaling.  This  would  be  the  same  effect  that  decoy  yields  for  DV.  The 
security  analysis  for  this  will  be  pursued  in  Phase  II. 


Figure  1:  Rate  vs  loss. 
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While  the  fundamental  rate  trade-offs  show  similar  trends  for  CV  and  DV,  CV  protocols  have 
lower  noise  (being  only  limited  by  local-oscillator  shot-noise  of  the  coherent-detection  receiver) 
and  can  access  a  higher  number  of  modes/second,  because  homodyne  or  heterodyne 
detection  can  potentially  have  a  much  higher  bandwidth  compared  to  single  photon  detection, 
at  comparable  detection  efficiencies.  On  the  other  hand,  error-correction  codes  are  better 
developed  for  small-alphabet  DV  protocols.  Because  DV  protocols  have  small  discrete  signaling 
constellations,  modulation  is  simpler,  as  compared  to  CV  protocols  (where  symbols  must  be 
chosen  from  a  Gaussian  distribution  or  a  densely-quantized  version  thereof).  Thus,  an 
extension  of  the  CV-BPSK  protocol  that  only  uses  a  few  modulation  constellation  points,  while 
achieving  the  R~r|  rate-loss  scaling,  will  ease  on  the  aforesaid  hardest  obstacle  to  CV 
implementations. 

Guha  et  al  have  invented  a  direct  communication  protocol,  that  is  quantum  secure  to  a  passive 
eavesdropper  (same  benchmark  of  security  as  the  Shapiro  two-way  protocol),  but  requires  only 
a  simple  one-way  binary-phase  laser-light  signaling,  near-LO-shot-noise-limited  homodyne 
detection,  and  a  reverse  authenticated  public  classical  channel  (which  may  be  an  RF  link  for 
instance).  The  bits/mode  performance  of  this  protocol  is  several  orders  of  magnitude  better  than 
the  Shapiro  protocol,  which  needs  entangled  states.  The  bits/mode  performance  achieved  by 
our  protocol  adheres  the  quantum-limited  rate-loss  scaling  (R~q),  and  is  only  factor  of  2  to  3 
below  it  for  reasonable  assumptions  on  sources  and  homodyne  detection. 

The  team  investigated  various  atmospheric  nonidealities,  including  the  effects  of  (1) 
atmospheric  absorption,  (2)  aerosols,  (3)  water  vapor,  (4)  turbulence-induced  amplitude  and 
phase  fluctuations,  on  loss,  and  that  of  the  blackbody  and  sky  radiance  on  detector  background 
counts.  The  main  objective  was  to  zero  in  on  a  few  candidate  windows  of  operation.  Taking  into 
account  the  trade-off  between  blackbody  radiance,  sky  radiance  and  atmospheric  transmission, 
the  three  candidate  wavelengths  we  identified  are  1.5  pm,  2.2  pm  and  4pm  (see  Figure  2).  Of 
2.2  pm  and  4  pm,  former  has  worse  (higher)  loss,  but  better  (lower)  noise. 


Figure  2:  Atmospheric  transmission  in  a  marine  environment  in  the  presence  of  water  vapor  and  aerosol 
(23km  visibility). 
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A  focal  point  in  Phase  I  was  our  investigation  of  the  performance  of  the  BB84  decoy  state 
protocol  at  these  three  wavelength  regions:  1.55  pm  ,  2  pm  and  4  pm.  Interestingly,  the  team 
found  that  in  the  presence  of  turbulence,  atmospheric  loss  and  sky/blackbody  radiance,  the 
three  wavelength  windows  exhibit  similar  performance  (Figure  3).  Thus,  the  choice  of 
wavelength  of  operation  really  depends  on  the  availability  of  optical  components  -  laser 
sources,  detectors,  modulators  -  at  the  various  wavelengths.  Given  the  easy  availability  of 
components  in  the  1.5  pm  telecom  band,  this  is  the  most  suitable  option. 
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Figure  3:  Performance  of  the  BB84-decoy  and  CV  Gaussian  modulation  protocols  across  wavelengths. 


The  Seakey  team  identified  the  1 .5  pm  FSO  communication  hardware  developed  by  Novasol 
Inc.  as  a  suitable  platform  from  for  implementation  of  a  QKD  testbed. 


Task  2.  Trade  study  that  includes  detailed  evaluation  of  atmospheric  turbulence,  and 
QKD  implementation  to  determine  optimum  operating  parameters  for  deployed  naval 
QKD  system.  Optimize  performance  across  choice  of  protocol,  code,  modulation, 
wavelength,  transmitter-receiver  technology. 

SoW:  Once  the  vulnerabilities  of  a  standard  free-space  QKD  system  operating  in  a  naval 
environment  have  been  identified,  Raytheon  BBN  will  design  a  QKD  system  complete  with 
technologies  to  mitigate  or  minimize  the  risks  incurred  by  operating  in  such  an  environment. 
Our  team  will  identify  the  QKD  protocol  (e.g.  BB84  with  decoy  states,  CV-QKD),  the  encoding 
and  modulation  format  (e.g.  polarization,  phase,  time-bin,  spatial,  including  high-order 
modulations),  the  wavelength  of  operation  (e.g.  specific  wavelength  within  the  visible,  NIR  or 
MWIR  band),  and  the  optimal  optical  channel  parameters  (e.g.  beam  waist). 
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Phase  I  Results:  The  performance  of  CV  Gaussian  and  DV  BB84  protocols  with  variation  in 
atmospheric  turbulence  is  shown  in  Figure  4,  which  plots  the  secret  key  rates  as  a  function  of 
channel  loss.  We  assume  the  fundamental  Gaussian  beam  is  modulated,  and  that  the 
apertures  are  circular.  For  this  plot,  all  the  “loss”  has  been  lumped  into  one  dB  figure  (the  x 
axis),  which  could  have  contributions  from  diffraction-limited  beam-spread,  atmospheric 
scattering  due  to  aerosols  and  water  vapor,  and  any  coupling  efficiency  loss  at  the  receiver — for 
instance  the  free-space  to  fiber  mode  coupling  efficiency  in  a  fiber-coupled  detector.  Realistic 
device  parameters  (as  listed  in  the  plot  legend)  were  chosen  for  both  CV  and  DV 
implementations.  As  expected,  the  key  rate  vs.  loss  degrades  with  increasing  Cn2.  Roughly 
speaking,  one  order  of  magnitude  increase  in  Cn2  results  in  the  key  rates  to  diminish  by  one 
order  of  magnitude. 


Figure  4:  Effect  of  turbulence  strength  on  secret  key  rate  of  the  CV  Gaussian  modulation  and  DV 
BB84-decoy  protocols. 

In  order  to  further  quantify  the  effects  of  various  marine  detriments  on  quantum  key  distribution 
rates,  our  team  built  a  numerical  model,  implemented  in  Matlab  (lnput_parameters.m).  This 
model  calculates  the  key  rates  of  the  BB84-decoy  and  CV-Gaussian  modulation  protocols  for  a 
specific  link  geometry  by  taking  into  account  user-defined  values  of  wavelength,  weather 
conditions,  visibility  (which  corresponds  to  aerosol  concentration  values),  and  elevation  above 
sea-level.  The  link  geometry  considered  accounts  for  transmit  and  receive-apertures  of  radius 
10  cm,  symbol  duration  of  5ns,  detector  efficiency  of  0.9  and  dark  count  probability  of  10'6. 
Figure  5  shows  a  screenshot  of  the  Matlab  GUI.  The  differing  weather/visibility  conditions  cause 
a  rapid  decline  in  secret  key  rate  of  both,  the  BB84-decoy  state  and  CV-Gaussian  modulation 
protocol.  As  a  best  case  scenario,  we  expect  key  rates  of  the  order  of  106  at  a  range  of  100m. 
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Figure  5:  Screenshot  of  Matlab  GUI 

We  further  examined  the  performance  improvement  by  using  multiple  spatial  modes.  We 
considered  a  soft-aperture  configuration  of  Hermite-Gaussian  and  Laguerre-Gaussian  modes, 
for  which  the  modal  transmissivities  are  exactly  solvable.  From  Figure  6,  it  is  seen  that  the 
advantage  of  using  multiple  spatial  modes  manifests  itself  at  shorter  ranges,  below  1  km.  At 
longer  ranges,  employing  multiple  spatial  modes  does  not  pose  any  advantage  due  to 
diffraction-induced  broadening. 
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Figure  6:  Secret  key  rates  as  a  function  of  range,  obtained  by  employing  multiple  spatial  modes 
(upper  figure).  Number  of  spatial  modes  as  a  function  of  range  (lower  figure). 


Section  C.  Publications 

•  M.  Takeoka,  S.  Guha,  M.  Wilde,  “Fundamental  rate-loss  tradeoff  for  optical  quantum  key 
distribution,”  Nature  Communications  5,  doi:10.1038/ncomms6235  (October  24,  2014). 

•  S.  Barzanjeh,  S.  Guha,  C.  Weedbrook,  D.  Vitali,  J.  Shapiro,  S.  Pirandola,  “Quantum 
Illumination  at  the  Microwave  Wavelengths,”  (to  be  published  in  PRL  in  March  2015), 
(January  31, 2015). 

To  be  published: 

•  S.  Guha,  M.  Takeoka,  H.  Krovi,  M.  Wilde,  C.  Lupo,  “Defeating  passive  eavesdropping 
using  a  laser-source  and  homodyne  detection”. 
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Section  D.  SEAKEY  Financial  Status 


Appendix  -  Summary  of  Phase  I  (PowerPoint  Slides) 
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Sea  Key  Goals 
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*  Evaluate,  adapt  and  optimize  free-spaceQKD 
technology  for  demonstration  and  deployment  on 
naval  platforms  such  as  surface  ships  linked  via 
free-space  channels  to  aircraft-based,  land-based 
or  other  ship-based  optical  transceivers 
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-  Identify  challenges  unique  to  marine  deployment,  rate- 
distance  tradeoffs  with  turbulence,  scattering,  absorption, 
background,  PAT  error  (e.g.,  due  to  platform  fluctuations) 
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-  Optimize  performance  across  choice  of  protocol,  code, 

"  modulation,  wavelength,  transmitter-receiver  technology 

-  Design  end-to-end  system  to  enable  key  exchange  over 
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L  simulation  of  reach,  rate  w/  naval  data,  implementation 
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-  Effect  of  turbulence  strength  to  key  rate  vs.  loss 

-  QKD  in  atmospheric  detriments:  wavelength  comparison 

-  Hardware  and  ancillary  technology  optimization 

-  Full  numerical  model  to  get  "parti als,:  to  key  rate  w.r.t. 
each  tunable  parameter  in  the  end-to-end  system 


QKD  over  single-modle  channel 
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Fundamental  tradeoff  between  QKD  key  rate 
and  input-output  power  transmittance,//  e  (0,  l) 


— Squash  *nle*igtofr>&nt  U6 

-  '  ■  Reverse  cohoMiri  nlccmalrai  LB 

-  ■  Wnnl  6(£W  r.octnt  V«y  rm 
-  BFSK  CV-CKD  ruwsu  rilfrnr— ftX  [xH»*sl 

Tnese  plots  as su  me  no 
excess  nose,  ideal  devices, 
no  tu  rbu  lence,  on  I y  sou  rce  of 
>n-deality  being  pure  bss 


Channel  kris  in  dl3  10  /??) 


-  UB  on  key  rate  for  any  QKD 
protocol, ipga(i  +  tf)/ 

k"'ol2  (  ^ 

meo^a,  Gu  ha,  Wilde,  arMv:  1310.0125  [quant-ph] 
submitted  to  Nature  Photonic*  -|  Match,  KU\ 

-  Ideal  BB84,  biased  bases, 
deterministic  SPS,.R(//)  «  lt 

Lo,  et  al.,  adC  ¥:quant-p-hi0011056  Journal  ot 
Crypto  og?  13,  Number 2;  133-16-5,  (2MS| 

-  Rev  coh-inf  LB.  achievable 
w  SPDC  +  displacement 
modulation,  l/loga(l  - 

Patron,  P1ra.ndola,  Llo^dt  Shapiro,  Phy£  Rev 
Lett  152,  215501  fiOOSj 


CV-QKD.  BPSK+ homodyne 


Lutnen  haul- group — PRA  73r  052316  2M6,  PflA 
76,  522313  2107  FRA  75,  012307  2C05 


Best-case  key  rate  for  point-to-point  QKD  must  decay  as R(q)  ^  ,  regardless 
of  QKD  protocol  or  power.  Optimal  transmit-powerdepends  on  loss.  SP-BBE4 
pretty  good.  Far  Held  LOS.-Dj  =  AtA^/iXL}2  <£  I  ,  transmittance^  k  Dj  oc  l/L2 
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Rate  vs.  loss 


& 


(R  oc  rr) 


[R  x  rt) 


10 


DV 

EES4- 

ro-cecff 


L-» 


cv 

EPSK 


ILo,  #t  all,  ar:rav:qua.nt-phrtHH1<Hfi,  Journal  oT 
Cl^ptMag}  13;  Num  trn  2,  134-164,  I2MSI 

1-ph  ideal  BB84 


J^ub  =  log 

I  —  TJ 

msoKa,  Guha,  Wild*. 
arXlv:  11115125,  to  appear  In 
Natural  Communication  |2G14| 

Patron,  Pinndala,  Lloyd, 
Shapiro,  PRL  102;  210401 


CV  (best 
achievable) 

=  hJK  f 


Zhao,  Hald,  R  gas 
Luti<.9nhaus,  P'RA  75, 
012307  <]2005| 


=  10" 

CV  \\=VJ  = 
Gafs&n  ;  i  DV 

-cc.  =:c-  i  *555fro&ar. 


20  40 

loss  in  dB 


6( 

10logi0{l/i7) 
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CV  and  DV  have 
same  optimal  rate¬ 
scaling,  R~n 

CV  binary  and  DV 
without  decoy  yield 
R~n2  scaling 

CV  with  few 
modulation  levels 
(not  full  Gaussian) 
should  retrieve 
optimal  (R~n}key 
rate  scaling  [same 
effect  that  decoy 
yields  for  DV],  will 
make  CV  a  LOT 
more  attractive 
security  pro  of  yet 
to  be  completed 
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-  Rate  vs.  loss 
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-  Direct  secure-communi cation  with  laser  light 

-  Effect  of  turbulence  strength  to  key  rate  vs.  loss 

-  QKD  in  atmospheric  detriments:  wavelength  comparison 

-  Hardware  and  ancillary  technology  optimization 

-  Full  numerical  model  to  get  "partials"  to  key  rate  w.r.t. 
each  tunable  parameter  in  the  end-to-end  system 
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CV  versus  DV 
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1.  Fundamental  rate-loss  scaling  [SAME] 

2.  Noise  [CV] 

-  CV  limited  by  LO  shot  noise 

-  CV  limited  by  single-mode  background 

3.  Signal  modes  per  second  (bits/mode-^  bits/sec)  [CV] 


-  Typically  driven  by  detection  bandwidth 

-  Homodyne/heterodyne  can  be  much  faster  Neither  of 

-  Wavelength  consideration  these  are 

4.  Error  correction  and  finite-key-length  [DV]  issues  against 

-  DV  better  studied  and  better  available  resource^>  cv-  EU 

5.  Ease  of  modulation  [DV]  made  goo™ 

-  Small  discrete  (phase/amp)  signaling  constellation  progresson4 

We  are 

6.  Repeater  [not  an  issue  for  this  project]  [DV]  considering  5 
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*  Summary  of  results 

-  Rate  vs.  loss 

-  CV  vs.  DV 

-  Direct  secure-communi cation  with  laser  light 

-  Effect  of  turbulence  strength  to  key  rate  vs.  loss 

-  QKD  in  atmospheric  detriments:  wavelength  comparison 

-  Hardware  and  ancillary  technology  optimization 

-  Full  numerical  model  to  get  ‘'partials”  to  key  rate  w.r.t. 
each  tunable  parameter  in  the  end-to-end  system 


High-rate  direct-secure  (LPI)  comnfcivttKwi 
using  DSSS  laser  light,  homodyneIBjTwhwk,g'L* 

*  Direct  quantum-secured  communication  immune  to 
passive  eavesdropping 

-  Does  NOT  need  entanglement  (as  in  Shapiro  protocol) 


X  |  =t  Or) 


r? 


tAlice 

Z  |  I  ±  ^ 

=  V  ©  £— 


rjai)  Eve 

t 


±  Homodyn^y'  =  X  ©  E 

I 

-W  =  Y  ©  V 


Bob 

G  u  ha ,  et  a  I . , 
unpublished 


-  Achieves  optimal  rate-loss  scaling:  within  factor  of  2-3  of 
best  possible  rate  (in  bits/mode) 

-  Shapiro  Ql  wproto col's  bits/mode  highly  suboptimal.  But 
very  high  modes/sec  due  to  the  naturally  broadband 
SPDC,  hence  resulting  in  a  higher  bits/second  than  above 


*  Can  we  get  around  this  shortcoming  of  our  protocol? 
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High-rate  direct-secure  comm.  Raytheon 
using  DSSS  laser  light,  homodyne>BjTwhwlqgiL* 

•  Direct  sequence  spread  spectrum  (DSSS)  CDMA  is  a 
way  to  tap  a  broad  BW 

•  Our  protocol  if  equipped  with  a  laser  transmitter,  AWG, 
and  a  broadband  homodyne  detection,  will  work  the 
same  way,  but  in  a  DSSS  mode  Will  achieve  MUCH 
higher  secret  bits/sec  than  Ql  protocol,  but  with  simpler 
hardware  (1  Mbps  over  10  dB  loss  should  be  possible) 

•  DSSS  is  a  also  popular  way  to  achieve  stealth 
(LPD/covert  communication),  and  anti-jamming 
capability 

•  We  recently  discovered  and  implemented  the  quantum 
limit  to  covert  communication  (sqrt(N)  bits  can  be  sent 
both  reliably  and  deniably  in  N  modes) 

BoulatA.  Bash,  Andrei  H.  Cheung  he,  Monika  Patel,  Jonathan  L.  Hahif,  Dennis  Goeekel,  Don  Tows  ley, 

and  SaikatGuha,  "Covert  Optical  Communication ",  in  review  with  Nature  (2tt1 4) 
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Atmospheric  LOS  turbulent 
channel:  spatial  modes 
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*  Focused  beam  $,(*■')  -  (  yr) ;  r  e  At,r'  e  a, 


3‘ 


-K-« 


Do  ±±  L 

1  (A  Lf 


k  =  2t/A 


-  Transmittance,  t/(£)  is  random  process  in  turbulence 


=  BfoM] ; 


L 


.  4JL2i  I/a 


[o^-Vx]  - 


\^HDr 


rh  D&)  =  i  .ug  k^C^L  xif* 

5hap  rc>  Ptiji  Rev  A  67  «2M$,  \%m\ 


-  Vacuum-prop,  channel  accommodates  infinite  orthogonal  spatial¬ 
mode  pairs.  Transmittances  add  up  to:  uy  =  y'^ 

-  In  turbulence,  at  any  time,  channel  can  accommodate  infinite 
mode  pairs.  Modes,  and  transmittances  are  random:  Dr  = 

*  Far  field:  one  usable  spatial  mode,  with  tji  as  Dj  1 

*  Nearfield:  ^  Dj  unity-transmittance  turbulent  modes 

*  Mean  transmittance  of  fundamental  mode  in  turbulence 
satisfies  the  following  bounds:  'te  <  £(fJii  ^  ,mnU>  Dj) 


_  Single-mode  turbulent  channel  Raytheon 

BBM  Technologies 


Mean  transmittance  of  a  Gaussian  FB 

Mean  transmittance  of  a  focused  Gaussian  beam 

-  First-order  estimate  of  SP-BB84  key  rate  in  turbulence,  using 
one  FB  (assuming  loss  and  turbulence  as  the  only 

*  Propagation  at  near-far  Held  boundary^ 
}  at  1 Os  of  km  range.  Multiple  spatial 
mod  e s  w  ill  h  e  Ip  o  n  ly  ma  rg  in  a  lly,  u  n  less 
la  rg  e  r  a  p  e  rtu  res  a  re  us  ed .  At  most  3-4 
modes  will  be  useful.  Orthogonal  spatial 
_  mod  e s  not  req  u  ired .  F  B s  s  u  fflcie  nt. 

'Dj  ;tt  \  —  riCIL 

at  A  =  Eti 

Dj  at-  A  =.  2  tti 

\  T-fJc;  at  CTjJ  =  10  Lr\  at  A  values  above 
(reasonable  assuTptfcn  fer  10-50  meers  awe 
sea  fesel  tirooeais  h  most  parts  o' fie  world: 

lrilomeler-B)  =  11)  -at  A  = 


.Wffi  1  GHz  pu  se 
proton  Cased 
CZf.  sflz  r  40C 
1  n -5  is  rg  5  cm  fac  ts 
1  u  ‘soever  aperture 
pots;  usilrg 
tu?1x  er  ce 


!0 


10 


Factors  that  deoreiEe  Hey-rate:  Enpbylrg  decoy-state  core  rent-state  E6S4  detector  eseess. 
noise  dart:  elilcfcs  laser  Intensity  noise  stfc-Lnity  coding  enclency  sneller  apertures  stronger  turbulence 
Factors  that  wil  boost  ley  irate:  using  nut  pie  (tew)  spat  bl  modes  both  poiarfcatlons.  larger 
apertures  v.esioer  (urtu  fence  snorter  waeergtr  :>«Ier20&e  more  hgrer-transmltarce  spatial  modes) 
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Secret  key  rate  (bits/s) 


Effect  of  turbulence  strength 


2  _ 


c 


10 


8  =  0,9  8  -  i 

rtd  =  0.9  Vd  =  o 
C  =  10  3  v  =  0. 

=  10  ,J  Nd  =  in 

CV  Gaussian  DV 
in  o  dulati  on  BB84~de 
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No 

scattering/absorpti 
on  inthese  plots 
CV  could  potential 
do  much  higher  rep 
rate  than  5  ns 
Lower  wave  lengths 
would 

accommodate 
multiple  spatial 
modes  (better  with 
CV)  but  worse 
turbulence 
performance  with 
fixed  channel 
geometry 


loss  iii  dB 
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Technical  update 

*  Summary  of  results 

-  Rate  vs.  loss 

-  CV  vs.  DV 

-  Direct  secure-communi cation  with  laser  light 

-  Effect  of  turbulence  strength  to  key  rate  vs.  loss 

-  QKD  in  atmospheric  detriments:  wavelength  comparison 

-  Hardware  and  ancillary  technology  optimization 

-  Full  numerical  model  to  get  "partials"  to  key  rate  w.r.t. 
each  tunable  parameter  in  the  end-to-end  system 


Sm X*y 


QKD  in  atmospheric  detriments 
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*  Narrow  down  wavelength  choices  using 
atmospheric  transmission  (absorption)  data 

*  Delineate  level  of  detriments  at  each  wavelength 

-  Loss:  contributions  from  (1)  atmospheric  absorption, 
(2)  aerosols,  (3)  water  vapor,  (4)  turbulence-induced 
amplitude  and  phase  fluctuations 

Michael  Jack,  Raytheon  RY5 

-  Noise:  contributions  from  thermal  bbn 

*  Additional  contributions  from  detector 

-  Loss:  detection  efficiency  (DE);  Noise:  dark  clicks  (Nd) 

-  Different  detector  type  (direct,  homodyne)  for  different 
protocol,  and  wavelength  of  operation,  affects  DE,  Nd 
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Loss  (1):  Atmospheric  absorption 
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Typical  curve  of  atmospheric 
transmission  in  clearweather 
[detector  response 
considered  Mndependent) 


H.  Man-3 r  and  S  Arnon,  AppL  C stc-E 
p.  am  ^2KHH| 


*  Why  not  operate  deeper  in  the  infrared? 

-  Quantum  cascade  lasers  (QCL)  capable  of  producing 
weak  coherent  states  between  3.5  pm  and  67  pm 

-  Problem  is  availability  of  single-photon  sensitive  detector 

Tam  porao,  ZWndsn,  Tanzllll. 

-  SFG-up-convert  to  visible  and  use  Si-APD  GlsBn- ■Afl,'9n.  Gkmannmi,  Faist 

r  d«  w^d,  ac,  no.  1-2  {2m\ 

■  SNSPD  [W-Si  or  other  material)  responsive  in  MWIR? 
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$ Mftjw 

Modeling  of  absorption 
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Mid-latitude  Summer,  elevation  10m, 
Ranges  of  1  km  and  7  km,  no  aerosol 


MrtKlUlh  SlHUNf 


Tropical  atmosphere,  elevation  10m, 
Ranges  1, 3, 7,  10  km,  no  aerosol 


*  RVS  used  MODTRAN  (Ontar  IP cModWi;n5,2)  software  to  generate  atmospheric 
transmission  estimates  for  MLS  (mid-latitude  summer)  and  tropical  environments, 
for  a  range  of  distances  atelevation  10m  above  sea  level 

*  The  absorption  bands  due  to  various  elements  in  the  atmosphere  [C02l  PhQ, 
ozone  etc.)  result  in  varying  transmission  levels 

*  Candidate  good  wavelengths,  purelyfrom  transmission  data:  1.55pm,  2.2pm,  4pm 

is  Hieraiu 


Loss  (2):  Aerosol  contribution 


ftayrtfietHi 
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*  Marine  aerosol 
concentrations  are  in  the 
range  of  100  to  300  cnr3 

*  The  size  distribution  usually 
contains  3  modes: 

-  nucleicD<100  nm 

-  accumulation 
1QQ<D<600  nm 

-  coarse  (D>600  nm) 


MalcriaJ 

JtefrO  | 

EmCnS 

duu 

W3 

I.S2  | 

S*UH 

] 

Walci  Kiluabk-T  walfr 

D.iH 

1.17 

$j6xio* 

T 

S «  <nlr  *■  motet  ("agedT  aeroMilj 

tili 

L.1& 

6.9*  ]0,j 

s 

SCi  Jilt  1  iv-H-TT  (tKW  tUOXt) 

2 

L.17  | 

6.5X3  01- 

4 

Sea  sail  ►  mulct  (nexr-uufKe) 

l  .It 

Tablr  J- AnTwl  m-aimal  HnnprciLian.  mem  radium  m:  J -d'r*.-lKc  indo  o'  Ihc  -t-niwiP- 
ANAM  amjwl  mwk?  fw  KAi  -  8G%.  Um  -  Uw  -  5  itvsk.  AMP  -  8*  h  ■  S  m, 

A  -  I.WSjua.  Phillip  Spfingla,  Joseph  Parana-,  NH.L,  2004 


*  Coarse  particles  makes  up  95%  mass 
but  only  5%  of  particle  number 

*  Sea-saltaerosol  concentrations  in 
marine  boundary  layer  (MBLJ—5-30  cm-3 

*  Aerosol  extinction  increasestowards 
sea  surface  due  to  greater  concentration 

*  Aerosoltransmission generally 
increases  with  wavelength  in  the  2-1 0 
pm  wave  length  range 


42  1  Z  1  i  4  ■*  T  *  9  IQ  M  iZ 


€.  Ai-naul  t’LcMmjhi  Fpvrtrn  Ttf  iJuTi-n-iil  Jur*  irffrnjgtM 
Aim*  h  b> vl>I  lift. 
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Raytheon 

Modeling  of  aerosol  contribution  80WT"^ 


Tropical  vs  Marine  Aero  Vis 


'.Vjviltripih  (um^ 

Elevation  10  m  from  sea  level,  range  =  7  km 

•  Marine  atmospheric  transmission  data  in  the  presence  of 
aerosol  for  a  few  candidate  visibilities  (10  km,  23  km  and 
50  km)  was  generated  from  the  MODTRAN  database 


Loss  (3):  Water  vapor 
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•  Combined  transmission  with  all  three  losses 

-  Mid-Latitude  Summer,  with  water  vapor,  maritime 
aerosol,  23  km  visibility,  7  km  range 


Midlatitude  Summer,  7  km  pathlength 
Total  tra  ns  w/ aero  and  water  vapor  continuum 


x  Good  transmission 


*  6 

WeiMfenp*  |unij 


windows  (1.55pm. 
2.2  ijm  and4pm ) 
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Loss,  noise  across  wavelengths 
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*  Note  that  of  2,2  |jm  and  4  ijm,  former  has  worse 
(higher)  loss,  but  better  (lower)  noise 


Radiances 

Sky  and  300  K  hl-adtbad^  id dianeeP 
1km  ae^asol  and  air  Eraritmiviirj-n 


k 


M 

P-H 

AT 

O.b  | 

I 

j 

AT 


9  7 
ftl. 
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Noise:  Background  and  Sky 
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*  Assumptions  for  background  count  calculations 

-  Spectral  filter  width:  1  pm 

-  Radius  of  receiver  aperture:  10  cm 

-  Field  of  view  of  receiver:  5  arc  second 


Wavelength 

Blackbody 

radiance 

(W  m'2  sra  |im'a) 

Sky  radiance 
(W  rrr2  sr1  |im*) 

Background 
counts/sec 
due  to 
blackbody 
radiance 

Background 
cowits/sec 
due  to  sky 
radiance 

1550  nm 

0.000000566 

1.5 

6.4-E-05 

169.6 

2200  nm 

0.000639 

0.2 

0.2 

32.3 

4000  nm 

0.7373 

0.009 

215.1 

2.6 
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SmMW 

BB84  across  wavelengths 


nayttieofi 

BBM  Technologies 


a 

"5 

in 

£ 


solid:  1.55^m  rt  -  tv  —  1Ucm 

dashed:  2.2[im  c*  =  w" 

dash-dotted:  4miti  T„  =-5  ng 

5  =  0.9  3  = 

t t,i  =  0  9  _  0.9 

f  “  L9^1e  V-  -  0.99 
t*  =  M>  b  ;V,(  -  10“' 
CV  Gaussian  DV 


10 


w 

o 

ryj 


Note:  To  be  fair  to 
10  PAT  capability, -one 
should  give  a 
diffraction-limited 
higher  aperture 
size  to  higher  A 

7  km  performance  points  include  losses 
absorption,  marine  aerosols,  and  '.vatervapo 
Bette  r  tra  nsmi  ssio  n  at  h  i  g  her  wavelengths 
a  I  mo  stco  mpensate  to  rwo  rse  mo  da  I 
^  tra  n  s  mitta  nee  i  n  tu  r bu  I  enc  e 


,fierosoi  absoqytkifi  and  scattering 
af feetstfie  transmission  coeff: 

overa  1 1  im  u  Iti  p  I  ier  to  kay  rata 

Sky. 'background  radiance  affect 
the  backg round  count  level: 

imposes  maximum  range 

We  considered  3  wavelengths, 
cased  on  transmission 


dulation  BBS4-decoy  background  tradeoff 


10 


■J  L  (in  kilometer b) 


10 


If  1-5  km  range  is  of  iiteiest  with  i=10cm  apertures,  or  for  a  30  km  range  vhnrtJh  if 
1=20011, 5-10  spat  a!  modes  worth  exploring,  when  at  1 .55pm.  At  2.2  pm  and 
4pm,  multimode  not  worthwhile  u n leas  40-5(teni  apertu res  a re  used. 
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Technical  update 

*  Summary  of  results 

-  Rate  vs.  loss 

-  CV  vs.  DV 

-  Direct  secure-communi cation  with  laser  light 

-  Effect  of  turbulence  strength  to  key  rate  vs.  loss 

-  QKD  in  atmospheric  detriments:  wavelength  comparison 

-  Hardware  and  ancillary  technology  optimization 

-  Full  numerical  model  to  get  "parti  als"  to  key  rate  w.r.t. 
each  tunable  parameter  in  the  end-to-end  system 


_ _  Hardware  and  ancillary  technology  navtiwon 

Cjb  8BN  Technologies 

optimization 

Identification  of  hardware  resources  for  stable  free-space 
transmission 

•  The  NovaSol  free-space  optical  communications 
interrogator  offers  position,  acquisition  and  tracking 
capability  at  1.5  pm  in  a  single  compact  box,  and 
has  been  tested  over  a  distance  of  50  km  in  a  naval 
environment.  This  operates  in  the  telecom-band 
region  of  1550  nm  and  can  be  modified  to  implement 
a  free-space  qkd  channel. 
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*  Summary  of  results 

-  Rate  vs.  loss 

-  CV  vs.  DV 

-  Direct  secure-communi cation  with  laser  light 

-  Effect  of  turbulence  strength  to  key  rate  vs.  loss 

-  QKD  in  atmospheric  detriments:  wavelength  comparison 

-  Hardware  and  ancillary  technology  optimization 

-  Full  numerical  model  to  get  "parti als"  to  key  rate  w.r.t. 
each  tunable  parameter  in  the  end-to-end  system 


Numerical  model 


Raytheon 

6BN  Technologies 


SeaJkey  Year  1  goafe: 

-  Identify  cna  enges  uniqu  e  to  raa r  ne  deployment,  rate-d  feta  nee 
tradeoffs wrtfi  turbulence,  scatter  ng,  absorptbn,  background,  RAT 
error  {e.g ..due  to  platform  fluctuations}. 

Optim  ee  performance  across  choice  of  protocol,  code,  mod  u  latbn , 
wavelength ,  transm  Itter-ieeeiver  technobgy. 


In  oider  to  address  tnese  goals,  we  have  fault  a  nu  merical  model  wtitidh 
generates  the  secret  key  rate  (m  farts,«see}  for  the  following  QKD 
protocols.: 

-  DV  BBM  with  aser  decoy 
CV  with  Gaussian  modulation 

This  numerical  model  Incorporates  following  user-defined  values : 
Wavelength 

Geog  raph  leal  zone  {eg .  Tropical,  Hld-lattude  etc.} 

Weather  eond  itbns  {eg .  clear,  cloudy,  raining} 

VTtsiifcflity  (co  nrespom  ding  to  a  sped  fie  ae  rosol  con  trifau  ton  } 

£  iev  at  on  a  cove  sea- level 

Tnese  values  are  drawn  from  the  MODTRAN  data  case,  wh  icti  is  an 
extensively  va  dated  atmospheric  radiative  transfer  model,  accessed 
fay  theOntar  Period  Win  fa.  2  software. 


INPUT 

QKD  protocol 
Gondftions  such  as 
geograpfi&al  location, 
range  ere. 


1  r 


Get  -p a rtials"  {sensitivity to  a 
parameteratan  operating  point) 


The  model  draws  the  transmission 
coefficientfromthe  MOOT  RAN  database. 


OUTPUT 

The  expected  secret  key 
.rare  fly?  fair ^mode)  as  drawn 
from  expected  values  of 
transmission.  background 
etc.  for  the  chosen 
conditions. 
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Instructions: 

-  Download  model.zip. 

-  At  the  command  prompt  in  Matlab,  type  Jlnput_parameters'.  This  will  run  the  m- 
code  which  produces  a  GUI  for  the  userto  enter  choice  of  wavelength,  weather 
condition  etc. 

-  The  file  !lnput_parameters‘  modelsthe  case  of  a  single  spatial  mode  and  takes 
turbulence  into  account.  A  second  m-file1“llnput_parameters_multimodes1 
models  multiple  spatial  modesin  a  soft-aperture  configuration,  and  does  not 
include  the  effects  of  turbulence. 

-  The  following  link  geometry  is  modeled: 


For  DV 


E  ftcror  cnosto1  tor ndyne 

d40CC30Ti  -  10*, ’Hoti* 
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Screenshot  of  Matlab  GUI 
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___L_  Taking  multiple  spatial  modes 
into  account 


Rarvtfiewi 
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Free-space  propagation: 

-  Farfield,  £>,  - 

'  £Al)» 


<  i 


-  One  spatial  mode  ( r?  ~  £>y)  has  any  appreciable  transmittance 

-  Rate.  R  -  l/L1 
-  Nearfield,  Df  1 

-There  are  ai  *s  d}  unit-transmittance  spatial  modes, i.e.,?jm  ^  1 

-  Rate,  R  -  l/L2 

-  Secret  key  rate  =  /?(x/])  +  2 Rfoi)  +  ...  +  M  ft(VM) 

■  If  we  consider  Hermite-Gaussian  and  Laguerre-Gaussian  modes  in  a 
soft-aperture  configuration,  their  modal  transmissivities  are  given  by: 

f  L  +  2Df  -  +  WS  Y 

^  w,  j 

Here,  Df  =  (kr?r/4L)  (kr2ltfl L)  isthe  productof  the  transmitter-pupil  and 
receiver-pupil  Fresnel  numbers  in  a  soft-aperture  configuration.  There  are 
s  p  ati  a  I  m  o  d  e  s  with  tra  n  s  m  iti  vitty  r . 
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Taking  multiple  spatial  modes 
into  account ...  contd. 
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A  few  examples 


The  adjacent  figure  depicts  the 
dependence  of  the  secret  key  rate  of 
the  BB 84-decoy  state  protocol  when 
multiple  spatial  modes  are  taken  into 
account.  As  expected,  this  is  more 
pronounced  atshorterranges. 

The  calculation  uses  a  series  of 
increasing  mode  numbers  until  the 
co  n  d  itio  n  w )/ ^  <  ,J- ]  i  s 

satisfied. 

The  rate  advantage  at  low  range  is  at 
the  cost  of  many  spatial  modes. 
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Solid  lines  UH-LaHudi  amie- 
dear  sip;  iDm  zJoae  sea-laei  aerosol 
vtetitmf  50  m.  ISHhii  rauefength 


Csshee!  lines:  HH-Ladude  Sunne- 
light  'itzzfe.  On  icw  ses-fte 
aerosol  SCtm  ISSOrm 

waeSff  gth 


ia 

a 

Gr 

€ 

f 

e 

D 


Du  E-E-S4  aftDC'j- 
GV-  Gjlss-  £r 

Modiiaticri 


-tzrge  (ffl) 


Transmission  quickly  drops  in  cloudy  and  rainy  weather  (from  =  0.96  to  >7=  0.5483 
and  =  0.16  for  light  and  heavy  rain  respectively,  leading  to  a  rapid  decrease  in  key 
rates.} 
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308d  limes:  MM-LattLrie  Summer, 
c  esT  st;  Cm  axue  seE-  a®  aerosol 
yfe-SW  ty  501m  Mrn  vrasergtr 


□ashed  iir^es:  Mld-L attune  Summei. 
dear  st;  Cm  eScm?  as-ftt  aerosol 
vlsMt;  SOm  400Crim  waeergir 


Taking  only  transmission  into  account,  the  rates  for  operation  at  2000nm  are  smaller 
than  those  at  4000nm.  Not  included  in  this  model  is  the  difference  in  detection 
efficiencies  and  background  count  rates  due  to  sky  radiance  and  blackbody  radiance. 


22^  A  few  examples 


contd. 
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SOW  mvz-  HH-Latfu*  Sumner, 
ceai  st;  Cm  sOae  as-fte  aerosol 
visJtfty  50  bn  ISSOth  ‘ffSK  ergtr 


Dashed  ires:  Mld-Lmide  Summer, 
deaf  st;  Cm  stxwe  =es-  ae  aerosol 
utsftlty  am.  ISSOrm  wasergr 


Range  :m;i 


50  km  visibility  is  relatively  high  for  coastal  environments,  and  this  only  occurs  on  very 
clear  days1.  If  the  visibility  drops  to  5  km,  the  key  rate  decreases  rapidly  as  shown. 


■VUsMtjj'  messu  remaps.  Song  esdended  pelt's  oer  ire  ccesn  surface.  J.  E  SMbUb  et  ai  Rroc  o'  SP  E  Vo!  5-5S1 


Raytheon  BBN  Technologies 

Distribution  Statement  A.  Approved  for  public  release ;  distribution  is  unlimited. 


Page  30 


gmmij 
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Including  multiple  spatial  modes(noturbulence): 


ScM  Wiio-LsatLce  Surma;  dear 
st;  Omaixwe  as  serose 

50  to  1550rn  ’Aaeergtr 


Dashed  ire  HH-LaHude  Suvvnef 
gft  ciEzJe.  On  ica^  see-  ae 
aerosol  ttsMy:  50tni  ISSCrm 
vraefength 


Tsrge  :mj 


Access  to  multiple  spatial  modes  at  low  ranges  results  in  an  increase  of  two  orders  of 
magnitude  in  key  rate. 


A  few  examples  ...  contd. 
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Including  multiple  spatial  modes  [no  turbulence}: 
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-  With  losses  from  aerosol  scattering,  water  vapor  continuum  absorption  and 
turbulence,  and  realistic  noise  levels  at  the  detectors,  the  secret  key  rates  of  the  CV- 
Gaussian  and  DV  BB-34  decoy  protocols  drop  rapidly  beyond  a  range  of  100-1000 
meters. 

-  For  the  lower  ranges,  the  ability  to  "turn  orf  multi-spatial-mode  operation,  and  choose 
up  to  a  certain  number  of  spatial  modes,  helps  in  increasing  key  rates  by  a  few  orders 
of  magnitude. 

-  It  may  be  more  realistic  to  assume  that  only  a  fraction  of  the  photons  not  collected  by 
Boh,  are  received  by  Eve.  We  will  in  corporate  this  correction  in  a  future  vers  ion  of  the 
numerical  model. 
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